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We study the excitation spectrum of strongly correlated lattice bosons for the Mott-insulating 
phase and for the superfiuid phase close to localization. Within a Schwinger-boson mean-field 
approach we find two gapped modes in the Mott insulator and the combination of a sound mode 
(Goldstone) and a gapped (Higgs) mode in the superfiuid. To make our findings comparable with 
experimental results, we calculate the dynamic structure factor as well as the linear response to 
the optical lattice modulation introduced by Stoferle et al. [Phys. Rev. Lett. 92, f 30403 (2004)]. 
We find that the puzzling finite frequency absorption observed in the superfiuid phase could be 
explained via the excitation of the gapped (Higgs) mode. We check the consistency of our results 
with an adapted /-sum-rule and propose an extension of the experimental technique by Stoferle et 
al. to further verify our findings. 

PACS numbers: 03.75.Kk, 39.25+k 



I. INTRODUCTION 

The prime example of a strongly correlated Bose sys- 
tem is 4 He — it exhibits normal, superfiuid, crystalline, 
and possibly even supcrsolid phases 1,2 ' '' and is attract- 
ing interest to this day. 1 Dilute cold bosonic atoms re- 
side generically in the weakly interacting limit but their 
tunability through quantum optical techniques allows for 
the realization of strongly correlated states. One way to 
achieve the strongly interacting limit is to reduce the 
kinetic energy by the application of an optical lattice, 
thereby effectively enhancing the effects of interactions. 
The broken translational symmetry then leads to new 
effects in the superfiuid phase not present in 4 He. The 
present paper is devoted to a study of the dynamical 
properties of lattice bosons within this strongly corre- 
lated regime. 

Jaksch et al ' pointed out that bosons in optical 
lattices are accurately described by the Bose-Hubbard 
Hamiltonian. Depending on the value of the nearest- 
neighbor hopping amplitude J, the on-site interaction U, 
and the chemical potential /i, the Bose-Hubbard model 
exhibits a superfiuid or an insulating (with uq £ N par- 
ticles per site) ground state, separated by a quantum 
phase transition.' The phase diagram has been investi- 
gated on a mean-field level, ' H using perturbation theory" 
and with numerical quantum Monte Carlo methods. 10 ' 11 
Apart from one dimension, the qualitative structure of 
the phase diagram is correctly described by mean-field 
calculations, cf. Fig. 1; however, fluctuations tend to 
shift the lobes to lower values of J/U in two dimensions. 11 

The dynamical properties of the Bose-Hubbard model 
were studied in both the Mott and the superfiuid phase. 
In the Mott phase, the excitation spectrum was stud- 
ied numerically, 9 ' 12 ' 1,3 on a mean-field level, 8 in a slave- 
particle approach, 14,15 and in strong coupling perturba- 



tion theory. 16 ' 1 ' All these approaches yield two gapped 
modes describing hole- and particle- type branches. For 
the weakly interacting superfiuid (J ^> U) the Gross- 
Pitaevskii ls ' 10 equation and Bogoliubov theory 20 pro- 
duce reliable results for the sound mode. However, for 
the superfiuid close to localization the character of the 
spectrum is still unclear. We will show below that this 
spectrum involves two distinct modes: (i) a sound mode 




FIG. I: The mean-field phase diagram of the Bose-Hubbard 
model involves disconnected incompressible Mott-insulating 
phases (grey shaded areas at small hopping J) where the 
density is pinned at integer values no, and a connected su- 
perfiuid phase in between and at larger values of J. The lines 
S^r mark the second order quantum phase transition sep- 
arating these phases. The commensurate filling in the Mott 
lobes is exported into the superfiuid along the lines &Hn (bent 
downward) where particle-hole symmetry is preserved. In our 
analysis we use a truncation scheme allowing us to discuss one 
lobe at a time. The black dots mark the positions in phase 
space where the spectra of Fig. 2 have been evaluated. The 
chemical potential difference Sfi is measured away from the 
lobe midpoint at J = 0. 
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characterized by a combined phase and density modu- 
lation, and (ii) a gapped mode describing exchange be- 
tween condensate and noncondensate at fixed overall den- 
sity 

The excitations of strongly correlated lattice bosons 
have been probed in two experiments. In the experiment 
by Greiner et al.' the whole lattice has been subjected 
to a global tilt. As this perturbation spoils the trans- 
lation invariance, it is difficult to infer bulk properties 
from their results; the observed resonant response has 
been discussed in Ref. 21. Second, in their recent ex- 
periment, Stoferle et al., 22 determine the energy absorp- 
tion due to a modulation of the lattice depth. They find 
a gapped continuum of excitations for small J/U and a 
broad feature in the supcrfluid. The appropriate response 
function has been calculated in the framework of the 
Gross-Pitaevskii equation 2 1 and for one-dimensional sys- 
tems (ID), both numerically 13 ' and analytically. 25 ' 26 A 
more powerful method to investigate this system is Bragg 
spectroscopy,"' described by the dynamic structure fac- 
tor 5(q, uj). As compared to the above lattice modula- 
tion technique, Bragg spectroscopy allows for a nonzero 
momentum transfer and hence reveals the full structure 
of the excitation spectrum. However, the nonzero mo- 
mentum transfer requires to apply two additional lasers 
at a finite angle to the system; so far, the limited opti- 
cal access has hindered an experiment using Bragg spec- 
troscopy on strongly correlated lattice bosons. On the 
theory side, the structure factor 5(q, u>) has been cal- 
culated in mean-field theory 28 and for ID systems us- 
ing numerical methods. 2f) ' !,) Below, we calculate both re- 
sponse functions within the Mott insulator as well as in 
the strongly correlated supcrfluid, making use of one uni- 
fied description. 

Our method is based on a technique previously intro- 
duced by Altman and Auerbach describing the particle- 
hole symmetric limit of the Bose- Hubbard model at large 
particle numbers no > 1. Here, we generalize this ap- 
proach to deal with the experimentally relevant regime 
of Mott insulators with one and two particles per lattice 
site. The method involves a truncation of the Hilbert 
space to three states per site and a spin-wave technique 
within a slave-boson language to describe fluctuations 
above a mean-field ground state. The truncation limits 
the validity of our results to a single Mott lobe and its 
surrounding supcrfluid environment. 

In the following Sec. II, we introduce the variational 
mean-field method outlined above. After proper trun- 
cation of the Hilbert space, we construct the mean-field 
phase diagram using a variational ground state. We pro- 
ceed with the derivation of an effective Hamiltonian de- 
scribing residual particle fluctuations using a method mo- 
tivated by the mapping to a spin-1 Hamiltonian. We 
diagonalize this effective Hamiltonian with the help of a 
generalized Bogoliubov transformation and find the spec- 
tra in both insulating and supcrfluid phases. Section III 
is devoted to the study of the response functions: we dis- 
cuss the dynamic structure factor (density-density corre- 



lations) in the Mott phase as well as in the superfluid and 
compare our findings with previously obtained theoreti- 
cal results. In addition to Bragg spectroscopy, we ana- 
lyze the lattice modulation technique which is described 
in terms of a dynamic modulation of the tunneling am- 
plitude J and calculate the corresponding response func- 
tion (hopping correlator) in both phases and for varying 
dimensionality of the excitation. Wc summarize and con- 
clude our work in Sec. IV. 



II. VARIATIONAL MEAN-FIELD 
A. Method 

The Bose-Hubbard Hamiltonian in a notation suitable 
for our approximation scheme takes the form 



#bh = -J^a\a 3 ■ + — ^Snj - Sfi^m, (1) 



where a] is the bosonic creation operator for a Wannicr 
state at site i, rij = a\ai is the number operator and 
Srii = m — no measures deviations of the particle number 
from a mean filling no. The chemical potential <5/i is 
measured from the middle of the lobe (cf. Fig. 1). 

Our goal is the determination of the dynamical proper- 
ties of the Bose-Hubbard model in the limit of strong in- 
teractions. In particular, we are interested in finding the 
excitation spectra and eigenstates in the Mott-insulating 
as well as in the superfluid phase nearby. Note that 
weakly interacting theories such as the Gross-Pitaevskii 
equation or the Bogoliubov theory cannot capture the 
physics close to localization. On the other hand, strong 
coupling pcrturbativc approaches are often incapable to 
correctly describe the broken U(l)-symmetry phase. 1 '' 
Altman and Auerbach !1 introduced a Hilbert space trun- 
cation method for large filling uq 1, i.e., the particle- 
hole symmetric case, which we extend to low fillings 
where particle-hole symmetry is broken. The basic idea 
is to truncate the bosonic Fock space to only three local 
states. In this truncated space, we first find a variational 
(mean-field) ground state and then derive an effective 
Hamiltonian H e fi for the excitations above this ground 
state. 

The truncation to three local states with particle num- 
bers no and no ± 1 is motivated by the strong suppression 
of particle number fluctuations in, and close to, the Mott 
phase; its validity is discussed in Sec. II C, below. We 
introduce bosonic operators that create 'particles' in the 
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retained three states 



*i,»|vac) = 



K T ) 



t\n +l 



VK + 1) 



vac) 



*-i,i|vac) 



(«l) r 



n : 



(2) 



where |vac) denotes the state with no particles present. 
The original bosonic operators ai can be expressed in 
terms of the t^j-operators (a = —1, 0, 1), 



(3) 



The Hilbcrt space spanned by the i Qi i-opcrators is too 
large and the physical subspace is obtained by imposing 
the constraint 



(4) 



a=-l 



The possibility to map the above truncated bosonic 
problem to a spin-1 Hamiltonian (sec Appendix A) moti- 
vates a strategy inspired by the spin-wave theory above a 
ferromagnetic or antiferromagnetic ground state: 52 start- 
ing out by expressing the spin operators S± and S z 
via Schwinger bosons a| B , bt B : S+ = as B 6 SB , S- = 
6s B a S B, and S z = (at B a SB - b\ B b SB )/2, the constraint 
at B a SB + bt B b SB = 2S is used in the ordered phase to 
go over to Holstcin-Primakoff bosons bl B — > b BP and 
aU -> V2S- btipbnp, with subsequent expansion of the 
square root in b B pb HP . In order to realize this program 
in the present situation, we first have to find the ground 
state playing the role of the ordered state in the spin 
problem. In a second step, we implement the holonomic 
constraint (4) via a procedure analogous to the change 
from Schwinger to Holstcin-Primakoff bosons. 

In order to find a proper ground state of the truncated 
problem we introduce the following variational operators 

b{i = 00809/2)4,. + sin(tf/2)[cos( X )iL + Mx)tl lti ], 
b\ , = - sin(tf/2)<, + cos(^/2)[cos( X )4 )i + sin( X )tl M ], 

(5) 



bli = - sm(x)t\i + cos(x)iLi 



ft 



where the Gutzwillcr-typc ground state shall be given by 



mx)) = IK> ac > 



(6) 



cf. (4). The relevant energy scale in the Mott phase 
is given by the interaction strength U. For the dis- 
cussion of the phase diagram all energies are rescaled 
and denoted with a bar, e.g., J = J/U. Combin- 
ing the Hamiltonian (1) and the ansatz (6) provides us 
with the variational energy per lattice site e yax (§,a) = 
a)\HBn\^{^, a))/N which is given by 



r (tf,cr) = sin(i9/2) 2 



1 



Sfi sm(2a) 



Jz 



sin 2 (tf) 



n + \J (no + l)^o cos(2cr) + -(1 + sin(2cr)) 



> (8) 



where z — 2d is the coordination number and N denotes 
the number of sites. In the Mott-insulating phase, par- 
ticle number fluctuations are absent within a mean-field 
approximation and thus $ — 0; the parameter a then 
drops out of the variational energy e vaY ($, c) and can be 
set to zero. In the supcrfluid case, i9 ^ turns out to be a 
convenient order parameter and hence a is eliminated via 
minimization of the variational energy (8) with respect 
to a, 



... 1 / 45/2 + Jz[cas(<d>) + 1] 

a[v) = - arctan — 

2 I 2Jz[cos(tf) + l]y/n{n+ 1) 



(9) 



This allows us to write e V ar(^;c) as a function of $ 
alone. Within a Ginzburg-Landau treatment of the phase 
transition, we reexpress $ in terms of the superfluid or- 
der parameter parameter ip = cr)\ai\^f , cr))/N = 
sin($)[^/no + 1 cos(7r/4— a) + ^/no sin(7r/4— cr)]/2 and ex- 
pand the variational energy (8) in ip, 



r ( t >, ( r( 1 ))) K a(J,W 2 + 



(10) 



The calculation of the coefficients a(J, Sp) and b(J, 8fi) is 
straightforward. The sign change of a( J, (5/2) marks the 
phase boundary and the roots of a( J, 5/2) = provide us 
with the known mean-field lobes 



Jz±dl- 2Jz(l + 2n ) + (Jz) 2 



(11) 



cf. Fig. 1. The tip of the lobes can be found by 
equating 5/2+ = 5/2~, providing the critical hopping 
l/(V"o" 



zJ c = l/(yno + 1 + v/no) . Due to particle-hole asym- 
metry, the line of integer density (a = 0) is bending down 
according to (cf. Fig. 1) 



The parameter ?9 controls the admixture of particle num- 
ber fluctuations in the ground state, whereas a deviation 
from integer filling is accounted for by a nonvanishing 
(7 = 7r/4 — x- As the transformation (5) is unitary, the 
b m ,i operators (m = 0, 1, 2) obey the constraint 



E = a, 

m=0 



(7) 



1 , 



<W = ~j[ Jz + W n o + 1 + V™oY 



(12) 



we refer to this line as the particle-hole symmetric line, 
which starts out from the tip of the lobe 5/2+ (J c ) as ex- 
pected. 

The above determination of the ground state (6) has 
provided us with the phase diagram of the Bosc Hubbard 
model and allows us to proceed with the second step of 
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our program, the implementation of the constraint (7) by 
going over to Holstcin-Primakoff-typc bosons; thereby, 
the operator i plays the role of the Schwinger boson 

ag B and the remaining operators b\ t , b 2i generate the 
excitations above this ground state, as does the operator 
&hp in the spin problem. We then eliminate one slave 
boson (6o,i) via the constraint 

bo,i = sjl - n^i - n 2 ,i, (13) 

where n TOj j = b^jpm^ (to = 1,2). Having chosen a 
good 'classical' ground state with potentially small fluc- 
tuations, we can expand the square root in (13) 

bo,% = \/l - ni.i - r%i ~ (1 - § n M _ i n 2,i); (14) 

the validity of the expansion will be discussed later (cf. 
Sec. II C). 

We express the Hamiltonian (1) in terms of the b 
bosons and eliminate the &o,i with (14). Collecting all 
terms up to quadratic order in the b m j provides the ef- 
fective Hamiltonian 

/ll,k f 12,1k. \ 

fl\ k fi k f 22 ;" k, (15) 

•lu.k ./12,k i/ll.k i/l2.k | 

f- 1 f- 1 q- 1 q- 1 J 

J 21, k J22,k y21,k i/22,k ' 

where 6k = (i>i,k, &2,ki &I _ k , ^2 -k) T anc ^ ^ denotes the 
first Brillouin zone. We have dropped a constant term 
to be discussed later. The coefficients of the normal g~* k 
and anomalous f~ 3 \ terms are given in Appendix C. The 
appearance of anomalous terms in the effective Hamilto- 
nian (15) is analogous to the situation where the corre- 
sponding spin problem is characterized by an antifcrro- 
magnetic ground state. These anomalous terms are re- 
moved via a Bogoliubov transformation which hybridizes 
creation and annihilation operators, thereby generating a 
new ground state carrying particle number fluctuations. 
This is in contrast to the ferromagnetic case, where the 
inclusion of quantum fluctuations does not impact on the 
classical ground state. The presence of the anomalous 
terms f~ B \ away from J = (cf. Appendix C), shows 
that the Mott-insulating state is 'nonclassical' and car- 
ries fluctuations for all finite values of J. 



B. Diagonalization procedures 

The diagonalization of the effective Hamiltonian (15) 
can be achieved via a (real) Bogoliubov transformation 

M/3 k = b k with D = M T H eS M (16) 

diagonal, where /3 k = (/3 K)k , f3 yM , fr x _ k , /3y,_ k ) T - In tlic 
Mott phase x = p (y = h) stands for particle and hole 
excitations respectively, whereas in the supcrfluid phase 
these indices stand for sound (x = s) and massive (y = rn) 



modes. Those bosonic commutation relations that are 
not automatically fulfilled are imposed by the additional 
condition 

MEM T = E, (17) 
where the matrix E is given by the outer product 

E= (& k 6+) T - (^) T (6 k ) T = diag(l, 1,-1,-1), 

reminiscent of the metric tensor in Minkowski space. The 
group 0(2, 2), given by all real 4x4-matrices M fulfilling 
(17), then shares many properties with the Lorentz group 
0(1,3), namely its decomposition in terms of 'boosts' 
(transformations of coordinates with different signs in the 
metric) and 'rotations' (transformations in a sector of the 
metric with equal signs). It turns out that this decompo- 
sition provides a useful strategy for the diagonalization 
of -f/eff in the Mott phase, where the symmetries of H c ff 
allow for an efficient determination of the corresponding 
rapidities and angles. On the other hand, in the super- 
fluid phase, these symmetries are absent and the diago- 
nalization of H e g is preferably done by mapping (16) and 
(17) to a non-Hcrmitian eigenvalue problem. 

(I) In the Mott state (cr = 0, -d — O), where no anoma- 
lous mixing terms f^ 2 \ = f 2 ~i\ = between the 61 )k 
and the 62. k bosons are present (see Appendix C), the 
parametrization of M in terms of boosts and rotations 
is suitable. To eliminate the anomalous terms and 
fmk one chooses a boost in the b m k -b m k plane; a subse- 
quent rotation in the &i, k -&2,k plane leads to the mean- 
field dispersions 

cp(h)(k) = Th(k)/2 + Sfj] +d>(k), (18) 

where w(k) = ^/U 2 - Ue Q (k)(4n + 2) + eg(k)/2 and 
eo(k) = 2J^ =1 cos(k • a;) is the bare band dispersion. 
Here, a; denote the vectors connecting nearest neighbors 
and we assume square and/or cubic symmetry, a = |a ( |. 
The dispersions shown in Figs. 2(a) and (c) characterize 
two modes, describing particle and hole type excitations, 
both with nonvanishing gaps A p (h) = e p(/i)(0), cf. Figs. 
2(e) and (f). In the Mott phase, the relation between the 
t a ,k and the b m ^ operators is trivial and the rotation in 
the 6i,k-fr2,k plane takes us back to the t a ^ operators. 
We therefore write the eigenstates in terms of the latter, 

/3j )k = A(k)4 k + S(k)t_ 1 ,_ k , (19) 
Pt M = -i4(k)tL lik -B(k)ti,_ k> (20) 

with A(k) = cosh(arctanh(/- 1 1 k /.g 1 ^ k )/2) and B(k) = 
sinh(arctanh(/ 1 " 1 1 k /g 1 " 1 1 k )/2), where we have used the re- 
lations g-^ = g- 2 ] k and f~] k = -f~%, which apply for 
the Mott phase (cf. Appendix C). 

(II) In the superfluid phase, the coefficients f~ 2 \ 
do not vanish and furthermore, /"^ 7^ — f^\. The 
presence of such terms renders a diagonalization via a 
parametrization of M as in the Mott state impractical. 
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FIG. 2: The panels (a)-(d) show the spectra for different points in the phase diagram (cf. Fig. 1) in the no = 1 lobe; all 
energies are measured in units of U. We display the results for the two-dimensional case with the dispersion along the direction 
k x . Panel (a) and (c) refer to the Mott phase, with the full line corresponding to the particle branch and the dashed line to 
the hole branch, (a) J = J c /3; black lines display dispersions on the line Sfi„ , gray lines correspond to <5/i = 0.2 U ; the change 
in chemical potential produces a shift in the spectra, cf. Eq. (18). (c) J = 0.9 J c and Su = 8u, l0 ; note the reduction in the gaps 
A p (h) as compared to (a). Panels (b) and (d) refer to the superfluid phase with sound- and massive modes given by dashed and 
dash-dotted lines, respectively, (b) J = 1.5 J c ; black and gray lines refer to the chemical potentials Su — 5fi no an d Su — 0.217. 
(d) J = 1.2 J c closer to the transition and with Su — 8fj, no . Panels (e) and (f) give the gap values A p , Ah, A m and the sound 
velocity v s for J £ [0, 2J C ] and Sfi = Su no (e) and Sfi = 0.2 U (f). The value J » 0.719 J c corresponds to the phase boundary 
for Su = 0.2 U (cf. Fig. 1). 



We therefore resort to the mapping onto a non-Hcrmitian 
eigenvalue problem. The constraint (17) written as 
M T = EM -1 £ and inserted into Eq. (16) yields 

M^EiTeffM = ED, (21) 

The problem of finding a matrix M £ 0(2,2) diagonal- 
izing H e g is now shifted to the problem of diagonalizing 
the non-Hermitian matrix EiT e g. The matrix M then 
is obtained from the eigenvectors {v} of Ei7 c ff via their 
proper normalization (with respect to E): let M be the 
matrix with columns {«}; then M £ 0(2, 2) is given by 

M = LM where L = diag(Zi, I2, h, h), 

with 




MEM T ) 

/ no 



(note that A-/EM T is diagonal, i.e., the eigenvectors {v} 
are automatically orthogonal with respect to the metric 
E). After diagonalization the effective Hamiltonian reads 

H eS = ^e s (k)/3t k /3 s>k + e m (k)/3^/3 m)k -C'^ £/ . (22) 
keif 

The optimization of the constant C 6 /^ leads to a renor- 
malization of 1? and a and is discussed below. The 
eigenvalues e s ( m )(k) can be calculated analytically (cf. 
Figs. 2(b) and (d)) and we find a sound (Goldstone) 
mode, which is linearly dispersing for k — ► with sound 
velocity v s = (9 fe e s (fc = 0)/ft and a massive (Higgs) mode 
with a gap A m = e m (0), cf. Figs. 2(e) and (f). Note the 



vanishing of the particle- and hole gaps and subsequent 
resurrection of the gap in the massive mode along the 
particle-hole symmetric line; for S/j, > 8fi no the hole gap 
transforms into the gap of the Higgs mode. The com- 
plete expressions for the dispersions of the sound and 
massive modes turn out to be lengthy and are given in 
Appendix C. 

To further characterize the excitations we consider 
coherent states of sound [|S s , q )] and massive modes 
[|#m, q )], respectively, 

|B. (Tn ), q )=e-l B .<»)-l a /V-^.^ t w-|0), (23) 

where |0) denotes the vacuum with respect to the 
/3 s ( m ) q opcrators (i.e., the new ground state) and 
SsM.q = l# s (m),ql cxp[i<5 s(m)iq ] arc complex numbers 
characterizing the coherent states. In Fig. 3 we present 
the expectation values of the density operator p, = 
(^s(m),q|ai a i|'8s(m),q) an( l the order parameter ipi = 
(# s (m), q |ai|S 5 ( m ), q ) for |B s ( m ), q | < 1 and q along x with 
q <C 7r/a, i.e., we only add a small amount of long wave- 
length excitations. The sound mode is given by a modu- 
lation of the phase accompanied by a modulation of the 
density, whereas the massive mode is given by a local con- 
version of condensate and noncondensatc. In an effective 
theory for the order parameter ipi, the sound mode cor- 
responds to the Goldstone mode, whereas the massive 
mode corresponds to the Higgs boson. The latter is only 
accessible for an effective theory which is second order in 
time. In Ref. 34, the author develops such an action for 
the superfluid to Mott-insulator phase transition. 

The eigenvectors can be calculated analytically as well, 
however, for our purpose the numerical solution of (21) 
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the broken-symmetry phase. In the Mott phase, C 5 j\j is 
not vanishing, but has its maximum again at d = a = 0, 
thus generating no renormalization. 
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FIG. 3: Expectation values for the modulus \ipi\ and the 
phase (j>i = arg^j) of the order parameter and the total den- 
sity pi for a coherent state of a massive mode (upper panel) 
and of a sound mode (lower panel). We have chosen param- 
eters no = 1, J = 1.2J C , Sfi = 5fi no guaranteeing particle- 
hole symmetry, resulting in a ground-state order parameter 
ipo ~ 0.47 marked by the faint gray line. Both modes refer to 
a q value along x and q <C n/a. The massive mode is given by 
a local conversion of condensate and noncondensate at fixed 
density, whereas the sound mode modulates the phase with 
a corresponding modulation of the density. In order to bring 
\ipi\ {4>i) in both panels in registry we have chosen a phase 



difference 5, 



= tt/2. 



is preferable. Hermiticity allows us to write the transfor- 
mation M in the form 



M(k) 



N(k) P(k) 
P(-k) iV(-k) 



(24) 



in addition, inversion symmetry renders the elements of 
the 2x2 matrices iV(k) and -P(k) independent of the sign 
of k and we can write 



4,k 



^i(k)/t,k + ^i2(k)< k + 
Pu(k)/? m ,_ k + P 12 (k)/Vk, (25) 
^2i(k)/3^ k + iV 22 (k)/3t k + 
P 2 i(k)/? m ,-k + P22(k)/3 s> _ k . (26) 



The matrices A^(k) and P(k) are calculated numerically 
and all quantities of interest, e.g., response functions, are 
given in terms of Nij (k) and Py (k) . 
The constant term 



U J,U 



E 



+ . 9 - 1 k -( es (k)+6 m (k))/2 (27) 



needs further discussion. The shift can be interpreted as 
a fluctuation-induced reduction of the ground-state en- 
ergy. In order to arrive at a self-consistent description, 
the parameters $ and a have to be determined by the new 
condition that the energy shift C 6 j\j is maximal. This 
corresponds to finding a new density and superfluid den- 
sity for a given set of parameters J, U, 5/x. Only if $ and 
a are renormalizcd with respect to their mean-field val- 
ues do we obtain a gapless sound mode, as demanded for 



C. Validity and expectation values 

The above method involves two approximations. First, 
the truncation of the Hilbcrt space: the quality of the 
truncation is expected to be acceptable in the Mott- 
insulating state and in the superfluid phase nearby, and 
turns bad in the weakly interacting limit. The qual- 
ity of the approximation can be checked through com- 
parison with experimental and numerical results, e.g., 
through testing the local number fluctuations. The lat- 
ter have been measured 55 and found to be suppressed 
due to strong interactions, at least in the vicinity of the 
Mott phase. The data is consistent with subpoissonian 
number statistics as predicted in Ref. 36 and supports 
the validity of the truncation to three local states. 51 The 
technique allows for systematic improvement by the in- 
clusion of further local states. 37 

Second, while the expansion of the constraint (14) is 
not a priori valid, we can check its quality via the (a 
posteriori) calculation of the ground state expectation 
values (6^ i b m J: we have found values in the range from 
0.18 (J — > 00) to 0.21 (J = J c ); fluctuations are largest 
at the phase transition and go to zero in the Mott phase 
at J — > for no = 1 which justify our expansion in (14). 

The calculation of matrix elements involves both states 
and the physical operator in question: so far, we have cal- 
culated the spectrum within the truncated Hilbert space, 
providing us with eigenvalues and eigenvectors of H c g . In 
addition, we need to express the operators in the eigen- 
basis generated by the (3 operators. The specific step 
of replacing 'Schwingcr bosons' by 'Holstein-Primakoff ' 
type operators involves the elimination of all &o,i opera- 
tors. Depending on the physical quantity under consid- 
eration, its expression through the (3 operators may or 
may not involve an expansion of the square root (14); 
in particular, the operators involving only bn, 62.1, and 
no.i = b\ x b 01 (e.g., the density operator in the Mott- 
insulating phase) can be transformed without requiring 
such an expansion. Otherwise, an additional imprecision 
has to be accepted due to the square root expansion of 
the constraint. 



III. RESPONSE FUNCTIONS 

While cold atoms in optical lattices excel in their tun- 
ability, they do suffer from a limited number of tools 
available for their characterization. In fact, so far only 
two experimental techniques are being used to determine 
the dynamical properties of cold atoms in an optical lat- 
tice, Bragg spectroscopy 2 ' and lattice modulation - ". In 
Bragg spectroscopy, two laser beams arc focused on the 
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system at an angle, leading to an inelastic two-photon 
scattering process. The system's response is described 
by the dynamic structure factor (density correlator) 

S(q,u) =^|(n|^ q |0)| 2 «5(^-^, lO ), (28) 

n 

where Spq is the density-fluctuation operator. The sum in 
(28) is running over all eigenstates \n) of the system with 
|0) denoting the ground state and fhv n Q the excitation 
energy associated with \n). This method provides anglc- 
resolvcd information on the system, a feature which, how- 
ever, turns out to be responsible for the method's limi- 
tation in actual experiments, as the optical access to the 
atom cloud is usually restricted. 

In the lattice modulation technique, recently intro- 
duced by Stoferle et ai, the depth of the optical lattice is 
modulated with a frequency u), introducing side bands in 
the laser forming the optical lattice. Within the frame- 
work of the Bose- Hubbard model (1), the lattice mod- 
ulation corresponds to a modulation in the hopping pa- 
rameter J; the determination of the energy transfer then 
boils down to a calculation of the hopping correlator 

§H = E \{n\T (x) \0)\ 2 5(hw - hj n0 ), (29) 



with T( x ) = } a l a j hopping operator (the in- 

dex x refers to a restriction of the modulation along one 
direction, here the x axis). The energy absorption rate 
then is proportional to toS(q, (j) and w5j^,"(w), respec- 
tively. In order to evaluate the response functions (28) 
and (29) we have to express the operators 5pq and 
in terms of the (3 bosons, involving a first transformation 
to b bosons and subsequent elimination of 6 ,i> cf- Sec. 
IIC. 

We first concentrate on Bragg spectroscopy in the Mott 
and superfluid phases. In the Mott phase, where the local 
density is pinned to an integer value in the ground state, 
we expect to excite a particle-hole continuum spread 
in energy as described by the bandwidth of these two- 
particle excitations. Surprisingly, we find pronounced 
peaks within this continuum which we can relate to the 
single-mode excitations &h(p) + e p(/i) s ) ■ I* 1 the superfluid 
phase one expects single-mode excitations, as breaking 
the U(l) symmetry is leading to collective excitations, 
and their weights will be determined. 

Second, we proceed with the calculation of the hopping 
correlator SY^(u>) in both the Mott-insulating and the 
superfluid phase. Again we find a continuum in the in- 
sulating phase. In the superfluid phase we do not expect 
to pump energy into the sound mode, as no momentum 
is transferred with this probe (up to a reciprocal lattice 
vector). A signal at finite energy will then give direct 
access to the massive mode. 



A. Structure factor in the Mott phase 

We make use of the eigenstates obtained within the 
variational mean-field approach to calculate the dynamic 
structure factor (28) in the Mott phase. In the truncated 
space the density-fluctuation operator 

*Pq = <a!«*»e- lq - r ' (30) 

i 

takes the form 

^P<± = *l,k*l,k+q — *-l,k*-l,k+q) 0*1) 

keif 

where we have used the exact constraint (7) to eliminate 
to,i (= bo,i in the Mott phase). Going over to j3 operators 
we obtain 

5 Pci =J2 A(k)B(k + q) [< k /C(k+ q) - 4,k/3;,_ ( k + q)] 

keA" 

and inserting this expression into the formula for the dy- 
namic structure factor yields 

S(q,w) = l [ — P(k,q)5[Hu;-e h (k)-e p (q-k)}, (32) 
z Jk v o 

with vo = (2n/a) d the volume of the Brillouin zone. The 
matrix element P(k, q) quantifies the coupling of each 
particle-hole excitation to the density perturbation and 
is given by 

1 + P(k,q)= (33) 
e (k)6o(q - k) + U 2 - U(2np + l)[e (k) + e (q - k)] 
4w(k)w(q-k) 

The dynamic structure factor (32) is closely related to 
the two-particle density of states (2DOS) given by 

r 

D(q,u)= — 5[hiu-e h (k)-e p (q-k)]. (34) 
Jk vq 

They share the same bandwidth as well as the character- 
istic total gap of the insulating state 

Atot = A p + A h = ^U 2 - 2JzU(2n + l) + (Jz) 2 . (35) 

Moreover, the van Hove singularities present in the 2DOS 
influence the response function, see Fig. 4(b). 

In our discussion of the results, 38 see Figs. 4(a) and 
(b), we concentrate on the two-dimensional situation 
where our mean-field analysis is sufficiently accurate. In 
the Mott phase, particle-hole excitations lead to a contin- 
uum, starting with a finite gap A to t at zero momentum 
q. The bandwidth decreases with increasing momentum 
transfer following the support of the 2DOS and reaches 
a minimum at the zone boundary, see Fig. 4(a). Figure 
4(b) shows the 2DOS D(q,oj) as well as S(q,u>) at q = 
(0.97r/a)e2;: The van Hove singularities in the 2DOS are 
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0.5 hu[U] 1 q x [1/a] 



FIG. 4: (a) Density plot of the dynamic structure factor 
5(q, w) in the Mott phase along q x for J — 0.9J C [point (c) 
in Fig. 1, energies are measured in units of U and no = 1, 
the chemical potential <5/i has no influence on the structure 
factor]. The spectrum is gapped, the spectrum's bandwidth 
decreases towards the zone edge, while at the same time devel- 
oping a pronounced structure allowing for the identification 
of single-mode excitations, (b) Cut along the uj direction at 
q x = 0.9-71 /a [dotted line in (a)]. The dynamic structure factor 
(solid line) exhibits marked peaks at the single-mode energies; 
the dotted line shows the 2DOS exhibiting van Hove singu- 
larities, (c) Single-mode weights of dynamic structure factor 
5(q, uj) in the superfluid phase at J = 1.2J C and 8[i = Sjj, no - 
The sound mode (Goldstone, solid line) exhausts all the avail- 
able spectral weight at long wavelengths, allowing the mas- 
sive mode (Higgs, dashed line) to gain weight only towards 
the zone boundary ([a.u.] denotes 'arbitrary units'). 



washed out in the response function which is dominated 
by the matrix element P(k, q). The latter generates the 
two pronounced peaks in 5*(q, uj) which we find (numer- 
ically) to be located at A p + €h{q x ) and Ah + e p (q x ). 
Hence, although Bragg spectroscopy generically excites 
a two-particle continuum, tracing these peaks allows for 
the extraction of the single-particle energies. These peaks 
are most prominent near the zone boundary and can be 
enhanced at small values of q x by dividing out a global 
modulation of the form [1 — cos{q x /a)\ (cf. Appendix B) 
from S"(q, uj). Furthermore, the peaks disappear deeper 
in the Mott phase where the excitations are more local- 
ized. We note that by expanding the integral (32) in J/U, 
we recover the results of the perturbative treatment (see 
Appendix B), hence the quadratic expansion of the con- 
straint (7) is consistent with second-order perturbation 
theory; however, the two results are comparable only for 
J < 10" 2 J c . 



B. Structure factor in the superfluid phase 

In the superfluid phase, the expression for the density 
fluctuation operator (31) written in terms of the b bosons 
contains bo.i operators, which have to be removed via an 
expansion of the square root (14), leading to terms linear 
in the operators b^ m k as well as higher-order terms. Going 
over to P operators we obtain the expression 

5 Pci = sm(tf/2){[P 21 (q) + iV 21 (q)]/4 )q 

[P 22 (q)+iV22(q)]/3l q }+0(/3 2 ). (36) 

Accounting only for terms up to first order, we ignore 
two- and multiparticle excitations and we cannot expect 
to fulfill the /-sum rule exactly, see the discussion in 
Sec. IIIC below. Inserting (36) into the expression for 
the dynamic structure factor yields 



^W(q, W ) = sin 2 (tf/2) 



S 2 (q)5(fuu-e s (q)) 



M 2 (q)6(tico-e m (q)) 



(37) 



this result reveals the two collective modes (22) charac- 
terizing the superfluid phase. Their weights are given 
by M (q) = P%\ (q) + N21 (q) for the massive mode and 
S(q) = P22 (q)+Af 2 2 (q) for the sound mode and are shown 
in Fig. 4(c). The sound mode is dominating the response 
at low momenta, with the massive mode acquiring weight 
only for higher momenta, where the sound mode satu- 
rates; in fact, numerical analysis confirms oc depen- 
dence of M 2 (q) at small q. 

Besides the single mode contribution, Spq also 
yields two-particle continua involving the excitations 
|m, k; s, q— k), |m, k; m, q— k) and |s, k; s, q — k). Their 
weight is about three orders of magnitude smaller than 
the single-particle contribution, however. 



C. Particle- number conservation 

The dynamic structure factor is constrained by several 
sum rules deriving from conservation laws. Gauge sym- 
metry and thus particle-number conservation is leading 
to the well known /-sum rule 



duj u)S(q, uj) 



Nq 2 
2m 



(38) 



which is modified in a one-band lattice description. 
The broken translational symmetry is leading to a 
nonquadratic dispersion, which can be characterized 
by a k-dependent effective mass tensor h 2 /m*j(]<.) = 
9fe i 9fc.eo(k). The /-sum rule 29 ' 39 adapted to the pres- 
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ence of a lattice then takes the form 

duwSfau) (39) 
^ {^o(k + q) + e (k - q) - 2e (k)}(0|a k a k |0) 

keA 

q-*o ^ q t qj ^ (Ola^aklO) 



Ef E 

( 



2 ^ m*Ak) 
tj keA v K ' 

Nq 2 



2m 



for 



m*(k) to 



Unlike in translation-invariant systems the structure of 
the ground state enters the /-sum rule via the nonuni- 
versal prefactor 



keA 



aj.dk |0) 



(k) 



(40) 



The expression (39) predicts oc q 2 behavior at small q, 
which is trivially fulfilled in the superfiuid phase (com- 
bine the weight S 2 oc q, see Fig. 4(c), with the linear 
dispersion of the sound mode; the oc q 4 dependence of 
M 2 does not contribute at small q) and can be easily 
verified in the Mott phase via expansion of the matrix 
element P(k, q), 

P(q,k)=pW(k)q 2 + 0(q 4 ). 

Unfortunately, our scheme does not allow for a precise 
calculation of the prefactor (40) and hence an exact self- 
consistency check (via particle-number conservation) of 
our result is not possible. 

The issue of number conservation has been raised in 
the work of van Oosten et al. 2 * . Their field-theoretic cal- 
culation of the structure factor did not reproduce the re- 
quired q 2 behavior, which then has been enforced through 
the use of Ward identities. However, it appears that the 
Green's function 8 G(iw„,k) used in this calculation al- 
ready violates number conservation in the Mott phase, 



iW ^G(i w „,k)# 
P Jk v o 



n . 



The application of Ward identities, although guarantee- 
ing number conservation, generates other defects in the 
structure factor, e.g., the appearance of linear terms in J 
spoiling the J — > — J symmetry present in bipartite lattice 
models. 



D. Lattice modulation in the Mott phase 

The lattice-depth modulation is a particle-number con- 
serving probe and hence produces only particle- hole exci- 
tations. While the lattice modulation has been uniaxial 



so far, 22 here, we also discuss its extension to an isotropic 
modulation (we discuss the case of equal modulation 
amplitudes but allow for mutual phase-differences be- 
tween the various directions). Expanding the constraint 
to second (i.e., leading) order, we obtain the response 
function 40 

^kin/, ,\ _ ± / D kin/ 



5 (x)M=o / — ^(kJ^-CfcCkJ-epC-k)], (41) 
z Jk v o 



with 



P ( k ^(k)=2n (n + l)|5] (:c) (k)| 



U 



(42) 



The interference generated by the different lattice mod- 
ulations is encoded in the sum 



£(k) = 2^V* cos(k-a ; ); 



i=i 



for the uniaxial modulation this reduces to 



S x (k) = 2cos(k -sl x ). 



(43) 



(44) 



The relative phase <fii between the different lattice mod- 
ulations can lead to interesting interference effects, see 
below. 

In Figs. 5(b) and (c) the results for isotropic and uni- 
axial lattice modulation arc shown for 4>i — and d = 2. 
The bandwidth is determined by the 2DOS (34) at zero 
momentum transfer q. Energy is transferred to the sys- 
tem only at frequencies to above the gap A to t/^j offering 
a simple way to determine the gap value. A dramatic 
change is obtained when going from the uniaxial to the 
isotropic modulation: the cusp at Hlo = U disappears and 
is replaced by a zero in the absorption probability. As 
the matrix element P kln (k) is non-negative, the response 
(41) only disappears if P kln (k) = for all k- values on 
the line defined by the <5-function in the integral (41). 
At nonzero 4>i, a finite weight is assembled away from 
the points (±7r/2, ±7r/2), leading to a finite response at 
cu = U. However, for <f>i tt an appreciable suppression 
is still observable. 

The result (41) obtained here has to be compared with 
the one obtained by Iucci and coworkers:" in their per- 
turbative calculation the factor [i7/d;(k)] 2 does not show 
up. Instead, the Bogoliubov transformation used here is 
equivalent to a resummation of diagrams and leads to 
this factor generating the interesting structure in the re- 
sponse function (41). 

Comparing our result with the experiment of Stofcrle 
et al. 22 we have to consider the uniaxial case, i.e., S^ ln (uj). 
While Stofcrle et al. observe a broad two-peak structure 
with maxima around U and 2/7, the present accuracy of 
our calculation does not account for high-energy excita- 
tions residing around 2U. On the other hand, the current 
experimental resolution does not allow to trace the inter- 
esting structure on the scale of the bandwidth, see Figs. 
5(b) and (c). Thus, both theory and experiment have to 
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0.2 hu [U] 1 1.4 0.8 hu [U] 1 
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FIG. 5: (a) Hopping correlator S^{ui) in the superfluid 
phase for J = 1.2J C and no = 1 for a parabolic trap and 
commensurate filling in the center (in a homogeneous sys- 
tem the response is nonzero only at uj = A m ). The inset 
shows the gap A m of the massive mode as a function of the 
chemical potential with a minimum at commensurate filling 
5fi = 5fj, no . The triangles are the calculated gap values and 
the solid line is a fit used in the calculation of S^ in (u) in a 
trap, (b) SfS(x) in the Mott phase at J = J c /3. The response 
consists of a particle-hole continuum with a gap A to t- The 
striking difference between the uniaxial [S* ln (u;), solid line] 
and the isotropic [S kln (cj), dashed line] situation is due to in- 
terference effects, (c) The same as in (d) for J = 0.9J C , where 
the gap Atot is larger and the available bandwidth smaller. 



be developed further in order to allow for a precise com- 
parison. Furthermore, it is worth mentioning that the 
current experiment may probes non-linear response. 41 



E. Lattice modulation in the superfluid phase 

In the superfluid phase (with an order parameter \ip\ > 
0), the expansion of the constraint again provides terms 
linear in the (3 operators and we obtain single mode peaks 
in the hopping correlator S^ ln at Huj = and hu = A m 
due to sound and massive excitations, respectively. En- 
ergy absorption is due to the excitation of the massive 
mode only and its weight W ma s S oc \ip\ 2 is given in Ap- 
pendix C. In addition, we also find a two-mode contin- 
uum (dominated by processes |s, k; m, — k)) with a weight 
suppressed by more than three orders of magnitude in the 
entire parameter range. 

In a trap, the sharp peak at A m will be smeared due 
to the inhomogencous density distribution. [Note that 
in the Mott phase, the response is not changed by the 
trap as both excitation energies e/j(k) + e p (q — k) (18) 
and the matrix elements P(k, q) and P^!"(k) are inde- 
pendent of Sfj] . The onset of absorption is determined by 
the minimal gap A™ n which occurs on the particle-hole 



symmetric line S/j, = 5/i no and the shape of the absorption 
profile depends on the distribution of effective chemical 
potentials /U e ff( r ) = Krap(f) m the trap. Here, we an- 
alyze the shape for a quadratic trap with /x e ff(0) = 5/U„ 
for a superfluid close to the Mott phase, J = 1.2J C . The 
dependence of the gap A m on the chemical potential <5/x 
can be calculated, see the inset of Fig. 5(a), and a con- 
venient fit is given by A m — A™ m = w(5n — 5fi na ) 2 with 
w ~ 2.06C/" 1 . The hopping correlator then draws its 
weight at frequency uj from rings matching the local gap 
energy A m (r), 



1 r R 

= ~m / dr27rrW mass 8[hu J ~A nl (r)} 

for fi^>A™ in . (45) 



1 



The resulting tail then resembles the broad absorption 
profile observed in the experiment of Stoferle et al. 22. 
Note that the precise shape depends on the actual density 
distribution in the trap; in particular, the divergence at 
A™ is removed when fx e e(0) ^ Sfj. n<3 . Another potential 
source of broadening is the finite lifetime of the massive 
mode due to the decay into two phonons as considered by 
Altman and Auerbach in Ref. 31. However, in two and 
more dimensions the effect of the trap dominates over the 
lifetime broadening. 



IV. SUMMARY AND CONCLUSIONS 

We have generalized the truncation scheme, intro- 
duced by Altman and Auerbach 1 ' to deal with the 
Bose-Hubbard model in the particle-hole symmetric limit 
(large particle numbers no S> 1), to the experimentally 
relevant situation of small filling numbers n of order 
unity. The determination of excitations consisting of 
small fluctuations about a variational ground state is in- 
spired by the Holstein-Primakoff description of quantum- 
spin systems and corresponds to the determination of 
spin-wave excitations above an antiferromagnetic ground 
state. 4 ' 5 We have determined the mean-field phase dia- 
gram as well as the spectra and eigenstates in the Mott- 
insulator and superfluid phases. These results then have 
been applied to the calculation of two response functions, 
the structure factor (density correlator) describing Bragg 
spectroscopy and the hopping correlator describing the 
lattice-modulation spectroscopy. 

A mean-field variational ansatz provides us with the 
usual phase diagram valid in dimensions d > 2. The in- 
clusion of two additional levels allows us to account for 
particle and hole-type excitations in the Mott phase and 
we determine, using a Bogoliubov transformation, the 
dispersions of the two gapped modes; at the same time, 
the Bogoliubov transformation introduces a new ground 
state carrying particle-hole fluctuations. In the super- 
fluid phase, we find a sound (Goldstone) and a massive 
(Higgs) mode and we determine the characteristic veloc- 
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ity and gap parameters. The massive mode describes a 
local counterflow of condensate and normal densities; this 
mode is absent in a Gross-Pitaevskii description where 
the dynamics involves a first-order time derivative, but 
is allowed in a Klein-Gordon type theory with a second- 
order dynamics. The presence of the latter is due to the 
underlying Mott physics providing 'particles' and 'anti- 
particles' (holes) and has been shown to be relevant in 
the Mott-insulator-superfluid transition in Ref. 21 (an 
additional first-order derivative is present away from the 
particle- hole symmetric line). A similar (massive) mode 
was found for the charge-density-wave compound NbSe2 
below the superconducting transition temperature. ' 
The fate of this gapped mode for U — > 0, where higher 
occupation numbers arc of importance, is currently under 
investigation. ! ' 

In the Mott phase, Bragg spectroscopy excites a 
particle-hole continuum and provides information on the 
gap and bandwidth of these two-particle excitations. To 
our surprise, we find that the structure factor unveils the 
single-particle excitation energies as well. In the super- 
fluid phase, collective modes (sound and massive) are ex- 
cited and visible as sharp peaks (to be smeared in a trap) ; 
the massive mode gains weight only at large momenta. 
The transition can be traced watching the appearance of 
a gap when crossing from the superfluid into the Mott- 
insulating phase. 

The lattice-modulation scheme 22 is presently the tool 
of choice to gain spectroscopic information on atomic 
matter in optical lattices. Our calculation of the hop- 
ping correlator providing the system's response reveals a 
two-particle continuum in the Mott phase which is sen- 
sitive to the details of the excitation scheme, uniaxial 
versus isotropic. The response in the superfluid is de- 
termined by the massive mode which appears as a sharp 
peak at finite frequencies. In order to better account for 
the experimental results we have extended our analysis 
to include the smearing due to the trap and find that 
the precise shape depends sensitively on the value of the 
chemical potential in the trap center. The experimental 
detection of such an energy absorption at finite frequency 
cannot be easily understood within a weakly interacting 
theory as described by the Gross-Pitaevskii equation. On 
the other hand, our strongly interacting theory provides a 
massive (Higgs) mode which naturally accounts for such 
a finite frequency absorption. Furthermore, a future ex- 
periment could address the question of how this Higgs 
mode disappears in the weekly interacting regime. 
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Kollath, F. Hassler, and A. Riiegg for discussions and the 
Swiss NSF for financial support via the NCCR MaNEP. 
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at ETH Zurich. Work in Innsbruck was supported by 
the Austrian Science Foundation, and the Institute for 
Quantum Information. 



APPENDIX A: SPIN HAMILTONIAN 

Defining the spin-1 operators , S~ and Sf in terms 
of the t a j operators via 

# = \/2(tI ji t ,i+4,i*-M) ( A1 ) 

ST = V2(4,i*i,i + *li,i*0,i) (A2) 

S' = 'i.,'i.< - tl lti t-i,i with (A3) 

[S+,Sr] = 2St, [Sf,St]=±St, (A4) 

we can write the Bose Hubbard Hamiltonian (1) in the 
truncated space as 



u 



(i,3) 



Jn Q £ 



^ ^ [s^ Sj S,i Sj~ -\~ Sj Sj~ ~\~ Sj Sj 



■ ttS*S+S-S* + S+S*S*S-)\ , (A5) 



where £ = \J (no + lj/no — 1 is a measure of the 'particle- 
hole symmetry- breaking'. In the work of Altman and 
Aucrbach, !1 £ was set to zero. 



APPENDIX B: PERTURBATION THEORY 

For a perturbative treatment of the dynamic structure 
factor one starts from the "pure Mott state |0) where all 
sites are occupied by exactly hq particles. A consistent 
expansion of (28) in J/U is obtained by an admixture of 
virtual particle hole pairs in the ground state 



i°> (1) — E i*> 



(i,j\aja m \b) 
U 



+ |6>, (Bl) 



where \i,j) denotes a state with no + 1 particles at i and 
no — 1 particles at j. The excited states in (28) are to 
lowest order given by the states 



with the energies 

e(k p , k h ) = U - (n + l)eo(kp) - n Q e (k h ). (B3) 

Inserting these perturbative states and energies into (28) 
leads to the expression 

SW{q,w)=N(j^ n (n + l)x (B4) 

[ * [d 7(q _ k) _ dl (k)] 2 5(hoj - e(q - k, k)), 
Jk v o 
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where 7(k) = 1/z ^ ;=1 exp(ik ■ a;). Both sides of the 
/-sum rule can be calculated independently and they co- 
incide with a value given by 



roc j2 

/ dujujS^(q,u) = 4— Ndn a (n + 1) 
Jo U 



i - ^7(q) 



which shows again the quadratic dependence on J and 
the vanishing with q 2 for small q. 



APPENDIX C: EFFECTIVE HAMILTONIAN 



1. Second order expansion 



Replacing all operators 6q.i m (1) expressed in the b- 
bosons and collecting all terms second order in bi t i, b%,i 
is leading to the effective Hamiltonian 



J 



keK 



2 i& I,k&l-k ■ 



(Cl) 



where we extracted a factor 2 Jz corresponding to the non-interacting band- width from the definitions of the coefficients 
given by 



9l2, k 

J ii, k 

r 1 

J 22, k 
J 12, k 



1 



1 



■ sin(2cr) cos(tf) + — cos(tf) + -[1 - cos 2 (tf)]{n + vM« + 1) cos(2cr) + -[sin(2<r) + 1]} 

£4 rj Z Q-tJ Z 



1, 



1 



I i[cos 2 W + l]{n+ ±[sin(2o-) + 1]} - ^vMn+l)[l - cos 2 (tf)] cos(2a) 



+ -jy-[3 - cos(tf)] sin(2cr) + -^-[cos(tf) + 1] + hi - cos 2 (tf))]{n + y/n(n + l)cos(2a) + hsm(2a) + 1]} 

"?{ n + ^ 1_sin(2(T)] } [c0S(,?) + 11 ' 

-^y- cos(i?/2) cos(2cr) + cos(i?/2)i[l - cos(tf)][cos(2er) + 2y/n{n + 1) sin(2a)] 

cos(t?/2) |[cos(t?) + 1] cos(2cr) + 2Vn(n+T)[l - cos(t?)] sm(2cr) j , 
y {[1 - cos 2 (i?)]{n + ^[sin(2t7) + 1]} - y/n(n + l) [cos 2 (tf) + l] cos(2ct)J , 



7k 
4 
7k 



-j s/n(n + l)[cos(tf) + 1] cos(2cr), 



s(i?/2) |[1 - cos(tf)] cos(2cr) + Vn(rT+I)[cos(tf) + 1] sin(2cr)| , 



(C2) 



where we defined again 7(k) = l/z^] 1=1 exp(ik • a;). 



2. Dispersion in the superfluid phase 



3. Massive mode weight the hopping correlator 



In the superfluid phase the excitation energies are 
given by 



e s(m) (k) = J Z] J2 [A e (k) T v / ^(k) 2 -4B e (k)j , (C3) 

and the coefficients are defined as 

A(k) = (.g-J 2 + (g-J 2 - (f~\f - (f~ 2 \f + 

2( ff -y 2 -2(/-:j 2 , 

B e (k) = ((^ - f-] k )( g - 2 ] k - f-\) - (sri - /r^) 2 ) x 



The weight W ma8S of the delta peak at A m in S km is 
^iven by 



W mass = sin(tf) 2 



cos( 1 9/2) 2 (A^ 21 (0) + P 2 i(0))- 



- cos(^)[V^+ Vm+T\ 2 (N n (0) + Pn(0)) 
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